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Abstract 

The nuclear liquid-gas phase transition of the system in ideal thermal equilibrium is studied 
with antisymmetrized molecular dynamics. The time evolution of a many-nucleon system confined 
in a container is solved for a long time to get a microcanonical ensemble of a given energy and 
volume. The temperature and the pressure are extracted from this ensemble and the caloric curves 
are constructed. The present work is the first time that a microscopic dynamical model which 
describes nuclear multifragmentation reactions well is directly applied to get the nuclear caloric 
curve. The obtained constant pressure caloric curves clearly show the characteristic feature of the 
liquid-gas phase transition, namely negative heat capacity (backbending) , which is expected for 
the phase transition in finite systems. 

PACS numbers: 24.10.Lx, 02.70.Ns, 24.60.-k, 25.70.Pq 
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I. INTRODUCTION 



Phase transition is an interesting phenomenon which appears in various physical systems. 
In nuclear systems with the excitation energies of a few to ten MeV/nucleon, the existence 
of the liquid-gas phase transition has been speculated based on the resemblance between 
the equation of state (EOS) of homogeneous nuclear matter and the Van der Waals EOS. 
However, the confirmation of the phase transition for realistic nuclear systems requires more 
careful arguments. 

First of all, any nuclear system accessible by experiments consists of a finite number 
of nucleons. Because of this, it may be generally believed that the phase transition tends 
to be smeared out due to the finite size effect, which is a correct statement for canonical 
ensembles specified by temperature. However, by studying microcanonical ensembles with 
fixed energies, a quite prominent signal of the phase transition is expected in finite systems, 
that is backbending (or negative heat capacity) in caloric curves 

HQS 

. Many experimental 

and theoretical works have been devoted to search for such a signal in various physical 
systems 

QBBBBB0Q. 

Another complexity comes from the fact that the excited systems for the study of nu- 
clear liquid-gas phase transition are produced in dynamical processes of nuclear reactions in 
laboratories. In fact, it is considered that the liquid-gas phase transition is somehow related 
to the multifragmentation phenomenon which is observed in the various nuclear reactions, 
such as in the relatively low energy region where the neck or midrapidity component may 
be created in dissipative binary reactions 0, 0, 0], in higher energy central collisions 



where clusters are produced copiously in expanding system 



15j . and also in peripheral col- 



lisions where the excited projectile-like fragment breaks up into pieces Thus intensive 
studies have been done to search evidences of the phase transition in these experimental 
data and some of the works have reported that indications of the phase transition have 

been obtained b y, h 0, 0, Q- 

However, the conclusion has not come yet and 
much efforts are still required. The difficulty is due to the nontrivial dynamical effects con- 
tained in the experimental data. Even though it may be true that the part of the system 



is equ i 



.i. 



21 



ibrated in good approximation as is indicated by the success of statistical models 



22, 23], the ambiguities enter in the data analysis through the identification of the 



equilibrated thermal source. On the other hand, there are several models which can predict 
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the experimental observables without such ambiguities by simulating the dynamics of these 



reactions microscopically 



24, 



25, 



26 



27 



28l |29[ . By carefully studying the time evolution of 
dynamical reactions, it may be possible to discuss how well the thermalization is achieved 
in dynamical collisions and how the results of dynamical reactions are related to the ther- 
mal properties, such as phase transition. For this purpose, the dynamical models should 
correctly describe the thermal properties as well as dynamical reaction mechanisms. It can 
be verified by applying these models directly to the ideal systems in thermal equilibrium, 
though the attempts for this direction have not been done sufficiently. 

The aim of the present work is to demonstrate that the nuclear liquid-gas phase transi- 
tion in ideal thermal equilibrium can be described by antisymmetrized molecular dynam- 
ics (AMD) which is a microscopic dynamical model based on the degrees of freedom of 
interacting nucleons [3, 13 13] ■ We utilize the same AMD model that has been ap- 



jlied to nuclear collisions and successfully reproduced various aspects of experimental data 



27 



3, H H 0, 3 0, S Q. 

To achieve this aim, the equilibrated systems are constructed and their statistical proper- 
ties are studied in the following way. Firstly, we prepare the initial state putting A nucleons 
(N neutrons and Z protons) arbitrarily within a fixed radius r wa u, and solve the time evolu- 
tion of the system by AMD for a long time. In order to confine the system in the container 
of volume V = |vrr^ all , we introduce a certain reflection process at the container wall. Then 
we regard the state at each time as a sample of the statistical ensemble with energy E, 
volume V and particle number A = N + Z. It should be emphasized that this ensemble is 
a microcanonical ensemble with a fixed energy E because the initial energy E is conserved 
through the AMD time evolution. By extracting the statistical information (temperature T 
and pressure P) from the ensembles, we can construct the caloric curves. The existence of 
the liquid-gas phase transition can be checked by finding a part with negative heat capacity 
(backbending) in these caloric curves, which is the characteristic feature of first order phase 
transition for finite systems Q, 

One of the advantages o: 
tional statistical models l],|21|, |22j, is that the calculation is based on the nucleon-nucleon 
interaction, rather than the nuclear binding energies and level densities given externally. 
Another advantage is that the interactions among fragment nuclei and nucleons are natu- 
rally taken into account, while they are usually ignored in other statistical models by the 



the present study as a statistical model, compared to conven- 



3 



freeze-out assumption. 

There have been several discussions on the question whether molecular dynamics models 
can describe the statistical properties of nuclear systems for which quantum and fermionic 



features are essential 



40. 



4l|. It has been shown that molecular dynamics 



can be consistent with the quantum and fermionic caloric curve if the wave function is 



ill 
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y antisymmetrized and an appropriate quantum branching process is taken into account 
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There arc several works 13 



38 



39 



41 



42, 



43 



44, 



451 ] that studied the nuclear liquid- 



gas phase transition using molecular dynamics models and some of the works claimed that 
a clear signal of phase transition, namely plateau or backbending in the caloric curves, is 
obtained. However, the molecular dynamics models utilized in these works have not been 
successfully applied to the dynamical multifragmentation reactions. Another problem is that 
the effective interactions adopted in these works sometimes do not satisfy the saturation 
property of infinite nuclear matter and none of them succeeded in obtaining plateau or 
backbending with an appropriate effective interaction so far. Furthermore, these works have 
calculated the constant volume caloric curves or those with a harmonic oscillator confining 
potential but have not explicitly shown the constant pressure caloric curves for which plateau 
or backbending is expected most clearly. On the other hand, our present study uses the 
same framework of AMD (AMD/DS) that has already been utilized for nuclear collision 

nn n 

simulations [30, 33] and the Gogny force [46J as the effective interaction between nucleons, 
which satisfies the saturation property of infinite nuclear matter. We also draw the caloric 
curves at constant pressure. 

This paper is organized as follows. In Sec. |H1 we show features of the Gogny force 
when it is applied to infinite nuclear matter. In Sec. IIIIl the framework of the AMD is 
explained, which is used to calculate the time evolution of many-nucleon system to create a 
microcanonical ensemble. The reflection process to confine the system in a given volume is 
also explained. In Sec. IIV[ the method to extract the statistical information (temperature 
T and pressure P) from our calculations is described. Then in Sec. the obtained caloric 
curves with constant volume and also with constant pressure are shown. The dependence 
of the caloric curves on the theoretical ambiguities is also discussed. Section EU is devoted 
to a summary and future perspectives. 
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FIG. 1: The EOS of symmetric nuclear matter. Each full line shows the pressure as a function of 
the density for the given temperature (0, 5, 10, 15 and 20 MeV). This EOS is obtained by adopting 
the Gogny force and assuming uniform density in the Hartree-Fock framework. Below the dotted 
line, the liquid and gas phases coexist if nonuniform density is allowed. The dashed lines are the 
EOS in this coexistence region obtained by the Maxwell construction. 

II. EFFECTIVE INTERACTION AND THE PROPERTY OF INFINITE NU- 
CLEAR MATTER 

In the present study, we adopt the Gogny force as the effective nuclear interaction. 
The Gogny force is one of the most successful effective interactions to reproduce the ground 
state properties of nuclei in mean field theories. The Gogny force satisfies the saturation 
property of the nuclear matter with the incompressibility K = 228 MeV when it is applied 
to infinite nuclear matter with the uniform Hartree-Fock approximation. The saturation 
property is essential for the discussion of the similarities and the differences between finite 
and infinite nuclear systems. The equation of state (EOS) of nuclear matter at finite tem- 
peratures can be also obtained by the uniform Hartree-Fock approximation, which is shown 
in Fig. ^ The dotted line is obtained by the Maxwell construction. Below this line, the 
system prefers to split into two parts, namely liquid and gas phases with different densities, 
rather than the uniform phase, and the actual EOS are the dashed lines in this coexistence 
region. 

The caloric curves for the infinite nuclear matter can be also obtained by the Maxwell 
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FIG. 2: The caloric curves of infinite nuclear matter obtained by applying the Maxwell construction 
to the EOS (Fig. The line of E/A = (3/2)T is drawn for the comparison, (a) The constant 
volume caloric curves drawn for the average densities p = OAOpo, 0.30/30, 0.23po, 0.18/30, 0.15po, 
O.lOpo, 0.07/3o, 0.04po, 0.02po and 0.015/9o, where po = 0.17 fm -3 is the normal nuclear matter 
density, (b) The constant pressure caloric curves for the pressures P = 0.02, 0.03, 0.05, 0.07, 0.10, 
0.15, 0.20, 0.25, 0.30 and 0.40MeV/fm 3 . The meaning of the diamond marks is explained in Sec. 

ED 

construction. The obtained caloric curves with constant volume and with constant pressure 
are shown in Fig. Efa) and Fig. Efb), respectively. The existence of phase transition can 
be recognized in both caloric curves. However, only the constant pressure caloric curves 
show temperature plateau at coexistent region. The constant volume caloric curves show 
monotonic increase of temperature with the energy increase even in the coexistent region. 
Finding a temperature plateau in the constant pressure caloric curves is an easy way to 
identify first order phase transition in infinite systems. 

In the realistic nuclear systems, however, the number of nucleons is several hundred 
at most. We have to carefully treat these systems because the ideas for infinite systems 
sometimes can not be applied to such small systems. Nevertheless, Gross et al. [H 0, [J have 
pointed out that the essence of phase transition is the anomalous concavity of the entropy 
S(E) = InW(E), with W(E) being the number of microstates, and the effect of phase 
transition can be clearly observed even in finite systems by finding backbending (negative 
heat capacity) in the microcanonical caloric curves. 
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III. FRAMEWORK 



The micro canonical ensemble with energy E, volume V and particle number A = N + Z 
can be obtained by solving the time evolution of the many-nucleon system confined in a 
container and regarding the state at each time as a sample of the ensemble. For the time 
evolution calculation, we utilize AMD, which is a microscopic dynamical model based on the 
degrees of freedom of interacting nucleons AMD uses a fully antisymmetrized 

product of Gaussian wave packets, 



(n...r A |$(Z)) = det 



exp< —v \ Tj — 



Z; 



XchU) 



(1) 



where the complex variables Z = {Zf, i = 1, . . . , A} represent the centroids of the nucleon 
wave packets. The width parameter v is treated as a constant parameter and taken as 
v = 0.16 fm~ 2 so as to reasonably describe the ground state of light nuclei such as 16 0. The 
spin-isospin states (xa, — P T> P I; n t; an d n j) are also independent of time. Because of 
the antisymmetrization, an AMD wave function |$(Z)) contains many quantum features so 
that it is even utilized for the study of nuclear structures [471 . 

The time evolution of the centroids Z, which parameterize |$(Z)), is determined by a 
stochastic equation of motion 

d 



dt 



Zi = {Z u H} + AZ t , 



(2) 



where {Zj,7i} is the deterministic term derived from the time-dependent variational prin- 
ciple and AZj is the stochastic term, which is important to describe the appearance of 
various reaction channels in multifragmentation, for example. The stochastic term AZ, ; is 
also essential for the consistency with the quantum statistics (37, 
of the stochastic term AZj are considered; one is the two-nucleon collisions and the other 
is related to the unrestricted single-particle motion in the mean field and the localization of 



38, 



40 



4l| . Two origins 



30 



32 



single-particle wave functions when fragments are formed 

In the following subsections, we briefly explain these terms referring to Refs. 
In this paper, we introduce several new improvements which are necessary for the present 
application because the energy conservation and the motion of emitted nucleons should be 
treated very consistently. 
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A. AMD time evolution 



1. Deterministic part of the equation of motion 



The deterministic term is derived from the time-dependent variational principle and given 



by 



where 



d 2 



with a,r = x,y, z 



27 



28, 



dZldZ ]T 



ln($(Z)|$(Z)) 



(3) 



(4) 



301 ] . The derived equation contains the expectation value of the 



effective Hamiltonian T~C, which is given by 



H(Z) 



mz)\HMz)) 



fc-mod(Z) + V mo d(Z). 



(5) 



(Hz)mz)} 

Equation (jSJ) contains the modification terms )C mo( i(Z) and V mo d(Z) for the kinetic and 
potential energies, respectively, both of which are for the treatment of the gaseous nucleons. 
If a nucleon is located out of nuclei and there are almost no other nucleons around it, 
the wave function of such a nucleon is interpreted to have a sharp momentum distribution 
rather than the momentum distribution corresponding to the Gaussian wave packet in Eq. 



This change of interpretation has 



nucleon emission and fragmentation 



27 



3een necessary for the consistency of Q-values of 



and is very important for the definition of 



temperature in the present work. By this modification, the zero-point kinetic energies of 



gaseous nucleons are subtracted from the total kinetic energy 

3ft 2 z/ 



27, m 



by 



2M 



■A + T (A-Af F (Z)), 



(6) 



where the zero-point kinetic energies of isolated fragments are also subtracted. The function 
Mf{Z) stands for the number of isolated fragments and nucleons which are assumed to 
have the definite center-of-mass momenta without zero-point energies. Therefore we take 
T (A— Mf) as the physical zero-point kinetic energy instead of (3h 2 i>/2M)A. The functional 
form of Mf{Z) is given in Appendix^! The parameter To should be 3h 2 u/2M in principle, 
but it may be adjusted to fine-tune the binding energies of nuclei. With the Gogny force, 
the binding energies for a wide range of nuclear chart are reproduced reasonably well by 
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taking T = 9.20 MeV which is close to 3h 2 u/2M = 10.0 MeV |3(j,|48(. In Appendix |H the 
degree of isolation X, for a nucleon i is introduced. Using Xj, the total zero-point kinetic 
energy can be written as 

A 

To(A-Af F )=T J2(l-Zi) (7) 

i=i 

so that the physical zero-point kinetic energy associated to the nucleon % may be considered 
as T (l — Xj). The potential energy should be also modified following the same change of 
interpretation of gaseous nucleons, because the spatial distribution of each gaseous nucleon is 
wide correspondingly to the sharp momentum distribution. The details of the modification 
term V mo a(Z) are also given in Appendix EI The parameters of the modification terms are 
chosen in such a way that a nucleon is treated as a gaseous nucleon (with a sharp momentum 
width) if the number of other nucleons around it is less than about 2. The binding energies 
of small nuclei such as deuteron and triton are utilized to fix some parameters. 

Equation (jSJ) defines the total energy that is to be conserved by the deterministic term 
{Zj, 7i\. The same energy is conserved also by the stochastic term. However, we need to be 
careful in treating the term JC mo d{Z) because the force — ^{Z^A/f} originating from this 
term is not always physically reasonable (see Sec. 1111 A"4|) . 



2. Stochastic treatment of single particle motion 



The stochastic term AZj in Eq. consists of two contributions. The first one is the 
stochastic two-nucleon collision term [23, 0,01 which is not shown explicitly in the following 
formulation for simplicity. The second one is related to the change of the phase space 
distribution by the mean field propagation. We follow the formalism of AMD/DS given in 



Refs. 



30 



but also introduce new improvements which are crucial in the present study. 
We write the Wigner function of the nucleon k as the mean of the stochastic phase space 
distributions of deformed Gaussian shape, 



f k (x,t) = g(x;X k (t),S k (t)) 



g{x;X, S)w k (X,t) 



d 6 X 



7T° 



(8) 
(9) 
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with 

g(x;X, S) 



x exp 



1 6 i 



2 

a,6=l 



We have introduced the 6-dimensional phase space coordinates 

x = {x } a= i,... i6 = jv^r, • (11) 

At an initial time, f k is represented by a single Gaussian wave packet g(x; X k , S k ) with 
Shn h = \5„h and the wave packet centroid X k is identified with the physical coordinate 
27. 



28j,|30j; 



X = {X fca } a=1 ,..., 6 = {ReW fc , ImW fc } . (12) 

For a short time, the time evolution of g k (x, t) = g(x; X k (t), S k (t)) by the Vlasov equation 

Sg^ _ _dh_ dg^^dli dg^ 
8t dp dr dr dp ' 

is characterized by the time evolution of the first and the second moments of the distribution 

d 6 3 



5 5 f d 6 x 

— X ka (t) = — / x a g k {x,t) — , (14) 

5 5 f d 6 x 

^S kab (t) = — (x a - X ka (t))(x b - X kb (t))g k (x,t) — , (15) 



in which (5/5t)g k (x,t) is given by Eq. ()13j) . A special notation of the time derivative (S/St) 
is adopted for the mean field propagation. The mean distribution [Eq. JSJ)] does not change 
when a part of the time evolution of the shape S kab is converted into a stochastic Gaussian 
fluctuation AX ka (t) to the centroid X ka (t) as 

-X fca (t) = j t X ka (t) + AX ka (t) (16) 



AX ka (t) = 0, (17) 



AX ka (t)AX kb (t') = D kab (t)5(t - 0, (18) 

where D kab (t) denotes the strength and correlation of fluctuations. Correspondingly, the 
equation of motion for S kab (t) is given by 

d $ 

J t s kab{t) = y t S kab {t) - D kab (t). (19) 
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In order that the Vlasov equation is satisfied, the choice of D ka b(t) is arbitrary as far as the 
positive definiteness of D ka b(t) and S ka b(t) is guaranteed. 

In the original version of AMD/DS in Refs. jsillssj], D ka b was chosen as the component 
diffusing beyond the original width of the wave packet, that is, 

= jim ^ £ max (o, A c - 1) O ac O bc , (20) 

where O a b and A a are the diagonalizing orthogonal matrix and the eigenvalues of the sym- 
metric matrix S k (t) + 4rSk At. With this original choice, however, there is no accordance 
between the momentum spreading of the distribution Sk and the zero-point kinetic energy 
T (l — I k ) assumed in the conserved energy (Sec. IIII A This accordance is important for 
the full consistency of the energy conservation and the precise evaluation of the temperature 
in the present work. Therefore we utilize the arbitrariness of D and introduce another choice 
D = + by adding a term as explained below. 

Ignoring the effect of antisymmetrization for simplicity, the expectation value of the 
kinetic energy of an emitted nucleon k can be written as 

2h 2 u 

^hamii = -JJj-( X li + X lr> + X ke) + T o(l - Ik), (21) 

which is a part of the conserved energy 7i and includes the zero-point kinetic energy. On the 
other hand, if we use the phase space distribution g(x; X k , Sk), the kinetic energy is given 
by 

2h 2 i> 2 2 2 4 
A^distr = ~jitf~(X k 4 + X kFj + X k6 ) + -T Tr p Sk, (22) 

where (?>h 2 v /2M) has been replaced by To consistently with Eq. (jSJ and a special notation 
Tr p Sk = Sku + Sk55 + Skm has been used. If only is taken, (1 — Ik) usually decreases 
faster than ^Tr p S k , when the nucleon k is going out of a nucleus. Therefore, in order to 
keep the accordance of /Chamii and /Cdistr, we choose D k so as to satisfy the condition 

Tr p S k = ^(l-l k ), (23) 

when the right hand side is getting smaller than the left hand side. This requirement can 
be satisfied by scaling S ka b as 

S kab (t + At) = R k S' kab (24) 
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by a factor 

fl ^ mm { ! W- 1 }- (25) 

where S' kab = S ka b{t) + {jiS kab - D^ b )At. This corresponds to taking the choice D kab (t) = 
D ( k l(t) + D®{t) with 

D k l(t)= Km o -^(l-R k )S' kab . (26) 

3. Decoherence process 

The fluctuation AX has been introduced in Eq. (jlfij) by the condition that the single- 
particle dynamics of the mean field propagation would be reproduced. However, the mean 
field propagation is not sufficient to approximate the time evolution of many-body systems, 
at least, in the sense that the idempotency of the one-body density matrix (p 2 = p) is 
unphysically kept during the mean field propagation. After enough time has past in the 
multifragmentation reactions, for example, the reduced one-body density matrix would be 
rather represented by an ensemble of density matrices in each of which single-particle wave 
functions are localized in fragments. This means that the coherence of the single-particle 
state is lost at some time (quantum branching). The decoherence is due to the many-body 
correlations and therefore its time scale (coherence time r ) should be related to many-body 
effects in some way. 

In Refs. I^nl . [jia], the decoherence is assumed to take place for a nucleon when it is 
scattered by a two-nucleon collision, which is the effect beyond mean field. In the present 
calculation, however, we choose a different prescription to investigate the dependence on the 
decoherence process. 

For each nucleon k, if there are more than three other nucleons within the radius of 
2 fm (measured in ReW/y^/), a decoherence process is assumed to take place with the 
probability of 1/tq per unit time. When a decoherence process takes place, it affects all the 
nucleons % that are located within the radius of 2 fm including the nucleon k itself and each 
of the single particle wave packet of the nucleon % is replaced by a Gaussian wave packet 
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with the phase space distributions 



S 



iab 



1 
4 







(a = b = 1,2,3) 
(a = 6 = 4,5,6) 
(a ?b) 



(27) 



The momentum widths are replaced by |(1 — Xj) rather than the standard Gaussian width 
| in order to satisfy Eq. (J2BJ)- It should be noted that a single decoherence process affects 
several nucleons at once and therefore the rate of decoherence for a specific nucleon is 
approximately proportional to the number of the neighboring nucleons. We take r = 500 
fm/c for usual calculations in the present work, but we also check the dependence of the 
results on r in Sec. IV CI 



4- The stochastic equation of motion 



The stochastic equation of motion for the wave packet centroids Z is given by 



HQ 



If 



Z % = {Z h H} +XX(AZ 4 )gg + ( AZ *)!f2p) ; (28) 



k=l 



(AZ i )g = {Z il O' k + T l k } Ck , (29) 



(AZO^MZi.WV (30) 

The fluctuation term (AZj)j^. is obtained by converting the fluctuation AX ka of the 
physical coordinate W to that of the original AMD coordinates Z. This conversion is done 
3Q by Educing a. expectation vatue O t( Z, t) o f a st oc h a stl c one- body ope rator 
which generates the fluctuation as {Z,,Cfc} by the Poisson brackets [Eq. Q]. 0' k in Eq. 
(I30|) is different from O k by the Lagrange multiplier terms for the conservation of the three 



components of the center-of-mass coordinate and those of the total momentum [30j,|32j]. The 
term T {Zj,X fc } in (AZj)j^ is introduced in order to cancel the unphysical force — T {Zj,X fc } 
in the deterministic term {Zj,7i}. Such a force does not exist for the single-particle motion 
in the mean field one-body dynamics. The subtraction is introduced as a term in (AZj)j^ for 
a technical reason (related to the dssp term) so that the Hamiltonian 7i with the zero-point 
correction term T (l — Ik) is still the conserved energy. 
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The term (AZj)^ is the dissipation term to achieve the energy conservation, where 
(Zia-jH 1 ) is defined by 

^n') = Wc^ ]T ^r. (31) 



J- 



Vflct 



The coefficient for each k is determined so that the energy violation by the (AZ 
term is compensated by (AZj)j^ . Lagrange multiplier terms are included in the effective 
Hamiltonian H', assuming that (AZj)^ p conserves some global one-body quantities such as 
the center-of-mass coordinate, the total momentum, the total orbital angular momentum, 
and the tnonopole and q uad rU p ol e tnotnents in the coordinate and tneanentntn 8 paees QQ. 

The subscripts and mean that the contents of the brackets are calculated by 
limiting to the subsystem or Cfc is the cluster that includes the nucleon k, where the 
clusters are identified by the condition that two nucleons i and j belong to the same cluster 
if |Zj — Zj| < 1.75. Nfc stands for a neighborhood of nucleon k defined by 

N fc = {i; \Wi - W fc | < 2.5, % EC k ,i^ k, and Mi > 1}, (32) 

where Mj is the number of nucleons within the distance of 3 fm (measured in ReW/^) 
from nucleon i. These limitations ensure that the fluctuation X^. a for the nucleon k should 
affect only the nucleons within the interactive range of the nucleon k when the conservation 
laws are imposed. The condition Mj > 1 in Eq. is introduced to exclude gaseous 
nucleons (here we regard the nucleon i with Mj < 1 as gaseous) from so that the 
dynamics of the gaseous nucleon i faithfully follows the single-particle motion determined 
by the deterministic term {Zi,H} and the fluctuation term ^fe=i(^^)fkt- The choice of 
Cfc and N fc here is similar to Refs. 3(1 , but has been updated in order to carefully treat 
gaseous nucleons, which is required for the precise definition of temperature. 

We cancel the fluctuation for a nucleon k when the number of nucleons in is less than 
five to prevent small clusters from unphysically breaking. 



B. Reflection at the wall 



To obtain a microcanonical ensemble with a fixed volume, we need to keep nucleons inside 
the container with a given volume V = ^r^ all during the time evolution. We introduce a 
kind of reflection process for this purpose, when nucleons or fragments are going out of the 



14 



container. At each time step after the time evolution without any effect of the container 
wall, we judge whether each nucleon k is in the container (^|ReW fc | < r wall ). If an isolated 
nucleon k is located outside the container and its momentum P k = 2hy/ulmZk directs out- 
ward, we change the momentum direction P& = Pfc/|Pfc| into an inward direction P' k which 
satisfies R& • P' k < 0, where Rfc = ^ReZ fc . We randomly choose the direction P' k as in the 
case of the reflection by an irregular surface. The absolute value of the momentum |P' fc | is 
adjusted so as to conserve the total energy which is sometimes affected by antisymmetriza- 
tion. The total angular momentum of the system is not conserved because of the irregular 
reflection, which allows us to construct a microcanonical ensemble without the constraint of 
the total angular momentum. 

When a nucleon which belongs to a cluster, where we regard the nucleon i and j belong 
to the same cluster when |W, — Wj| < 0.8, is located outside the container, we apply a 
similar reflection procedure to the center of mass coordinate of the cluster. If one applied 
the reflection procedure to eacfl of the nucleons in the Ouster as is done in Ref, Q flS 
14 (1 14 1L 1441 . 145| , the cluster would be unphysically broken by the crash to the wall. 

When a nucleon k is reflected, the shape S k has to be reflected consistently with the 
change of the momentum direction Pfc — > P' k . We consider the reflection with respect to the 
plane that includes the point Rfc and is perpendicular to n*, = (P^ — Pfc)/|P' fc — P&|, and 
hence the component parallel to rifc is reversed. According to this operation, the coordinate 
or momentum vector in the intrinsic frame of the wave packet k is transformed into 

x' fe = Xfc - 2 (xfc • rifc) h k = Tfc x fe , (33) 
(Tfc)cTT = 5 aT - 2n ka n kT (a, r = x,y,z). (34) 

By defining the transformation matrix in the phase space as 

%=( Tt °), (35) 
the transformation of the shape Sk is given by 

S'kab = '^2(.?k)ac(Tk)bdSkcd- (36) 

cd 

In the actual calculation, for the purpose of numerical stability, we introduce a small 
delay time rdei ay for the response to the fluctuation AX ka j3, Q • This delayed fluctuation 
Sfe a (i) also has to be transformed as E' ka = J2b(1~k)ab'^kb when a nucleon k is reflected. 
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IV. CALCULATION OF TEMPERATURE AND PRESSURE 

In order to obtain caloric curves, we need to extract the temperature T of each micro- 
canonical ensemble {E, V, A}. We also need to know the pressure P of the ensemble in order 
to draw the constant pressure caloric curves. In this section, we explain how to obtain T 
and P in our calculation. 

The micro canonical temperature is defined by T" 1 = (dS/dE)\v t A- This quantity can 
be evaluated with a good precision by using the average kinetic energy of gaseous nucleons. 
When the total volume is V > A/p , where po — 0.17 fm~ 3 is the normal nuclear density, 
it is possible to find a partial system composed of gaseous nucleons to which classical ideal 
gas relations are applicable. Let us call this partial system G which will be defined below. 
Appendix [B] shows that the microcanonical temperature T of the ensemble {E, V, A} can 
be obtained from the average kinetic energy of G-subsystem as 



where JCq and Aq are the kinetic energy and the number of nucleons of G-subsystem, 
respectively, and the brackets ( ) {e,v,a,a g >o} denote the average value for the microcanonical 
ensemble with Aq > 0. Because of the consistency of Eq. (J2T|) and Eq. (J22J), A^g can be 
calculated as 



The first term is the contribution from the wave packet centroid. The accumulated delayed 
fluctuation rdei ay Sfe |22| should be included in this term because the corresponding shape 
shrinking has been already applied. The second term is the contribution from the momentum 
widths of Sk- 

According to Appendix El G-subsystem can be defined arbitrarily as far as the following 
two conditions are satisfied. The first condition is that the quantum effect should be negligi- 
ble for any nucleon in G-subsystem. The second condition is that the G-subsystem should be 
selected without using momentum variables, i.e., the configurations with different nucleon 
momenta should be equally taken into account if they have the same nucleon positions. 

In our actual calculation, G-subsystem has been chosen in the following way. We select 




(37) 




(38) 
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the nucleons k for which the density of the nucleons with the same spin-isospin 



P°= E ^] 3/2 e- 2 ( ReWfc - ReW ^ 2 , (39) 



is sufficiently low (p a < pc = (l/200)po) so that the antisymmetrization effect can be ne- 
glected. It is necessary to eliminate the nucleons belonging to clusters because the discrete 
level(s) of the internal degrees of freedom can only be treated quantum mechanically. There- 
fore, among the selected nucleons, we choose the nucleons k which do not have more than 
one other nucleon within the distance of tq = 3 fm (if we use defined in Sec. IIII A 21 we 
choose the nucleon k which satisfy the condition < 1). Because of our reflection proce- 
dure, it is possible that nucleons locate outside the container in short time interval. Those 
nucleons are excluded from the selected nucleons. The above selections are done based on 
spatial coordinates without any momentum selection. We define the system composed of 
these selected nucleons as G-subsystem. The results should be independent of the definition 
of G as far as the necessary conditions are respected, which will be checked in Sec. IV Gl bv 
changing the criteria in our definition of G. 

For the calculation of the pressure P, we adopt the common definition of pressure as the 
external force necessary to keep the volume. Namely the pressure is given by 

P = An/r ttl £ A P'N ( 4 °) 

OT «D TtaW reflections 

where the summation is taken over all the reflections at the container wall (see Sec. IIIIB|) 
which occurred during the time evolution of the total time r to tai, Ap is the momentum 
change at each reflection, and N is the normal vector. Factor two is from the fact that when 
a nucleon or a fragment hits the wall the total momentum of the rest of the system is also 
changed for the momentum conservation. 

V. THE STUDIED SYSTEM AND OBTAINED CALORIC CURVES 
A. The studied system 

In the present study, t he sy stem with (N, Z) = (18, 18) is considered, which is the same 



lesy, 

Q. 



system taken in Refs. |44l. |45| . There are several examples indicating that phase transition 



exists even in such small systems (4J, |5| 
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5 10 15 20 25 30 

E*/A [MeV] 

FIG. 3: The constant volume caloric curves for the A = 36 system obtained by AMD. The lines 
correspond to the container size r wa ii = 5, 5.5, 6, 6.5, 7, 8, 9 and 11 fm from the top. Statistical 
uncertainty is shown by error bars. The line of E/A = (3/2)T is drawn for comparison. 

Micro canonical ensembles are constructed with the radius of the container r wa n = 
5, 5.5, . . . , 15 fm, which correspond to the average densities A/(^-r^ an ) = 0.40po) 
0.30p , • • • , 0.015p , and the energy E*/A = 4, 6, . . . , 28 MeV, where E* stands for the 
excitation energy relative to the ground state of 36 Ar nucleus (E gA = — 8.9 A MeV). The 
minimum radius r wa n = 5 fm corresponds to a volume of the container V = 2.5A/p , and 
therefore we do not consider a compressed liquid nucleus in this study. 

The time evolution was calculated up to 20000 fm/c, which is much longer than a typical 
nuclear reaction time scale (~ 100 fm/c). The states for the first 5000 fm/c were discarded 
in order to remove the initial state dependence. Samplings were done at every 10 fm/c. 
Similar calculations were carried out 4 times independently to improve the statistics. 

B. Obtained caloric curves 

The obtained constant volume caloric curves with r wa u = 5, 5.5, 6, 6.5, 7, 8, 9 and 11 
fm are shown in Fig. |3J The caloric curves with r wa n = 13 and 15 fm are also constructed, 
but they are not shown in Fig El because the complete equilibration could not be achieved 
by our investigation time (20000 fm/c) in the low energy cases (E*/A = 4 ~ 8 MeV) 
[the results systematically differ if we start the calculation with a quite different initial 
state]. Otherwise, fairly smooth caloric curves are obtained. The negative heat capacity 
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FIG. 4: The constant pressure lines drawn on E-V plain. The circles indicate the points where 
the microcanonical ensembles are constructed. The lines correspond to the constant pressure lines 
with P = 0.02, 0.03, 0.05, 0.07, 0.10, 0.15, 0.20, 0.25, 0.30 and 0.40 MeV/fm 3 , which are obtained 
by performing the transformation between volume V and pressure P with the interpolation at each 
energy E. Statistical uncertainty is shown by error bars. 

(backbending) or plateau does not appear in the constant volume caloric curves, which is 
consistent with the result of Ref. 45] with another version of AMD. It was expected from 
the fact that the constant volume caloric curves for infinite nuclear matter do not show the 
plateau even in coexistent region (Fig. 12(a)). 

Although this result (Fig. EI) is somehow similar to the caloric curves for infinite matter 
(Fig. 12(a)), the transition from the liquid-gas phase coexistence to the pure gas phase is 
not clear. It is difficult to judge clearly whether the phase transition exist or not from the 
constant volume caloric curves. On the other hand, the constant pressure caloric curves 
which showed plateau in the case of infinite nuclear matter (Fig. 12(b)) are expected to show 
the signal of the liquid-gas phase transition most clearly. Therefore, the judgment should be 
done with the constant pressure caloric curve, not with the constant volume caloric curve. 

Constant pressure caloric curves can also be constructed from our results as follows. 
After the evaluation of all the ensembles {E, V, A}, we eventually know the temperature 
and the pressure at each lattice point on the E-V plain indicated by circles in Fig. 0J From 
these results, we can draw constant pressure lines on the E-V plain by performing the 
transformation between volume V and pressure P with the interpolation at each energy E. 
The obtained constant pressure lines are shown in Fig. 0] The relation between energy E 
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FIG. 5: The constant pressure caloric curves for the A = 36 system obtained by AMD. The lines 
correspond to the pressure P = 0.02, 0.03, 0.05, 0.07, 0.10, 0.15, 0.20, 0.25, 0.30 and 0.40 MeV/fm 3 
from the bottom. Statistical uncertainty is shown by error bars. The curves of E* j A = T 2 /(8 MeV) 
and E*/A = T 2 /(13MeV), and the line of E/A = (3/2)T are drawn for comparison. 

and temperature T along these constant pressure lines corresponds to the constant pressure 
caloric curves. These caloric curves should not be confused with the caloric curves for the 
constant pressure ensemble {E,P,A}. The obtained caloric curves are those calculated 
from the microcanonical ensemble {E,Vp, A}, where the volume Vp is chosen so that the 
ensemble gives a certain pressure P. Figure E] shows the constant pressure caloric curves 
with the pressure P = 0.02, 0.03, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30 and 0.40 MeV/fm 3 
obtained by the above procedure. In each caloric curve with P < 0.15 MeV/fm 3 , negative 
heat capacity is observed clearly, which is the signal of first order phase transition. For 
example at P = 0.05 MeV/fm 3 , the caloric curve start from E*/A = 8 MeV and the 
temperature decreases till the energy of 16 MeV even the system is heated up from 8 MeV 
to 16 MeV, and after the energy of 16 MeV, the temperature goes up with a slope 3/2. 
By the comparison with the infinite nuclear matter caloric curves (Fig. I2^b)), the caloric 
curves with P < 0.20 MeV/fm 3 in Fig. El can be interpreted as the caloric curves drawn 
from the liquid-gas phase coexistence to the pure gas phase. Figure El shows the fragment 
mass distributions for the ensembles along the P = 0.05 MeV/fm 3 line. Full lines are the 
distributions of the fragments identified by the condition that two nucleons % and j belong 
to the same cluster if |^Re(Wj — Wj)| < r c i ust = 2.5 fm. These distributions do not 
necessarily correspond to the fragments which would be emitted when the container wall is 
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FIG. 6: The fragment mass distributions along the P = 0.05 MeV/fm 3 line. Full lines are the 
distributions obtained with r c i ust = 2.5 fm. Dashed lines are the distributions obtained with 
?~ciust = 3.0 fm. 

removed. In fact, the distributions change when r c i ust is varied. The dashed lines in Fig. 
|B]show the fragment mass distributions obtained with r clust = 3.0 fm. Nevertheless, these 
distributions are helpful for the qualitative understanding of the change of the system with 
energy increase. When the energy is low (E* j A = 8 MeV), the distribution shows U-shape 
with two peaks so that the typical configuration at this energy is a large nucleus coexisting 
with a few gaseous nucleons. When the energy is increased [E* j A = 12 ~ 16 MeV), the peak 
at the large fragment becomes smaller and the distribution changes into shoulder-like and 
power-low-like shapes. Thus complex configurations with many intermediate and light mass 
fragments are typical at these energies, and the proportion of light fragments increases as 
the energy increases (from E* / A = 12 MeV to 16 MeV). When the energy is sufficiently high 
{E* j A = 20 MeV), the distribution changes into exponential shape, which can be interpreted 
that the nucleons are moving almost freely at this energy, although a few nucleons come 
close and make compounds with small probability. The change of the distribution is fully 
consistent with the interpretation of the caloric curve, that is, the system along the P = 0.05 
MeV/fm 3 line changes from liquid- gas phase coexistence at low energy to pure-gas phase at 
high energy. 

The caloric curves in Fig. El do not contain the pure liquid phase because of our choice of 
the container size r wa n > 5 fm. Nevertheless, we can roughly guess the pure liquid caloric 
curve from our result as follows. Diamond marks in the constant pressure caloric curves 
for infinite nuclear matter (Fig. Efb)) indicate the points corresponding to the density with 
r waii = 5 fm. These points are inside of the plateau region, but are very close to the left edge 
of the coexistence region. Furthermore, as is seen in Fig. Efb), the pure liquid caloric curves 
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FIG. 7: The temperatures obtained with different parameter sets (rQ and pq) for choosing G- 
subsystem. The results of the systems with r wa n = 5,9,15 fm and E*/A = 6,14,22 MeV are 
shown. The abscissa {Aq) is the average number of nucleons in G-subsystem and the parameters 
are chosen so that (^4g) becomes about half and about a quarter of the original choice [tq = 3 fm 
and pg = (l/200)po]- The open points show the dependence of tq and the solid points show the 
dependence of pq, which are connected by lines to guide eyes. Statistical uncertainty is similar to 
or smaller than the size of points. 

are almost independent of the pressure unless the liquid is compressed with higher pressure 
(P > 0.3 MeV/fm 3 ). We can guess the corresponding line by connecting the leftmost points 
of our AMD result and the line can be compared to the Fermi gas formula E* /A = aT 2 for 
excited nuclei. The pure liquid caloric curve should be slightly left to the connected line. 
Thus our AMD result seems to be consistent with the quantum and fermionic statistics with 
a reasonable level density parameter a -1 = 8 ~ 13 MeV. 

If we simply define the critical point as the point where the negative heat capacity disap- 
pears, the critical temperature and the critical pressure can be estimated as Tq ~ 12 MeV 
and Pq ~ 0.20 MeV/fm 3 , respectively. These values are lower than those of infinite nuclear 
matter {Tq ~ 16 MeV and Pq ~ 0.3 MeV/fm 3 ). The lowering of Tq seems to be reasonable 
because the existence of the surface for fragments reduces the interaction energy gain by 
being fragments (liquid phase). 
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C. Check of the theoretical ambiguities 



As is mentioned in Sec. IIVI the measured temperatures should be independent of the 
choice of G-subsystem as far as the conditions described in Sec. IIVI are satisfied. This can 
be checked by changing the criteria in the definition of G. In Sec. I1V[ G-subsystem has been 
selected as the subsystem of the nucleons which satisfy p a < Pq = (l/200)po an d do not 
have more than one other nucleon within the distance of vq = 3 fm. 

According to Appendix El the dependence on pq and tq should be weak if only exist. 
This is actually demonstrated in Fig. [7| which shows the temperatures obtained with different 
choices of parameters (tq and po) to select G-subsystem. The parameters are chosen so 
that the average number of nucleons in G-subsystem (Aq) becomes about half and about 
a quarter of the original choice. The open points show the dependence on tq and the solid 
points show the dependence on pq. The temperatures are almost independent of the choice 
of the parameters. Only small dependence on the parameters can be noticed when (Aq) < 1. 
This is probably because we measure the temperature by using only the states with Aq > 
(see Appendix |H}. Nevertheless, the dependence is very weak and seen only at low energy 
and at small volume so that the ambiguity associated with this issue would hardly affect 
the obtained results shown in Sec. IV Bl and would not change the conclusions. 

In the current AMD dynamical calculation, the coherence time To (Sec. IIII A3|) is some- 
how arbitrary. The above results (Figs. Old) are obtained with the choice r = 500fm/c. We 
should check the robustness of the obtained results against the arbitrariness of the coherence 
time. Thus we have done the calculations with r = 250 and 1000 fm/c. The obtained con- 
stant pressure caloric curves are shown in FigOJ^a) and Fig. Efb), respectively. Furthermore, 
we have done the calculation with another choice of the decoherence process adopted in Refs. 
(30! I33I ] and the obtained constant pressure caloric curves are shown in Fig. IBJc). By the 
comparison of the figures in Fig. |H] and Fig. OJ we notice some change of the leftmost point 
of each caloric curve, which can be interpreted as the change of the level density parameter 
a for the liquid caloric curve. This dependence limits the possible range of the coherence 
time r when the level density parameter a is known. In any case, negative heat capacity is 
observed in all the caloric curves for a wide range of r . 
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FIG. 8: The AMD results of the constant pressure caloric curves with P = 0.02, 0.03, 0.05, 0.10, 
0.15, 0.20, 0.30 and 0.40 MeV/fm 3 from the bottom when the coherence time is varied: (a) tq = 250 



fm/c, (b) ro = 1000 fm/c, (c) the same procedure as in Refs. 
VI. SUMMARY 
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In this paper, we have constructed the microcanonical caloric curves of a system with 36 
nucleons [(N, Z) = (18, 18)] using the microscopic time evolution of AMD. We applied the 
same AMD model that has been applied to nuclear reactions successfully. 

For the first time as a molecular dynamics calculation for nuclear systems, we have con- 
structed the constant pressure caloric curves and found that negative heat capacity appears 
there with an appropriate nuclear force (the Gogny force). Negative heat capacity (back- 
bending) is a specific character of first order phase transition in finite systems. Thus we 
confirm the existence of nuclear liquid-gas phase transition under the ideal thermal equi- 
librium condition with AMD. The obtained fragment mass distributions also support the 
existence of the phase transition. We have checked that the extracted temperature from the 
ensembles do not change when different criteria to choose the nucleons for the temperature 
measurement are used as long as necessary conditions are satisfied. We have also checked 
the obtained conclusions are insensitive to the theoretical ambiguity, i.e. the coherence time 
t . The critical temperature and the critical pressure are estimated as Tc ~ 12 MeV and 
P c ~ 0.20 MeV/fm 3 , but these values may slightly depend on the choice of ro. As for 
a future study, it is interesting to investigate the dependence of the backbending and the 
critical point on the size and the isospin composition of the system. 

We have also shown constant volume caloric curves, but the signal of the phase transition 
can not be seen clearly there. It suggests that it is important to identify the path of caloric 
curves drawn on E-V plain when those are constructed from experimental data to discuss 
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the existence of nuclear liquid-gas phase transition. 

As we have demonstrated in this paper, AMD can be applied to ideal thermal equilibrium 
situations as well as the dynamical processes of the reactions. Therefore, it may be possible 
to construct the method to extract the statistical information from experimental data by 
studying the reaction simulation results of AMD together with the calculations of ideal 
systems in thermal equilibrium. 

Acknowledgments 

This work was partly supported by High Energy Accelerator Research Organization 
(KEK) as a supercomputer project. 

APPENDIX A: THE MODIFICATION OF THE ENERGY FOR GASEOUS NU- 
CLEONS 

The fragment number Mf utilized in this paper is similar to the one in Ref. but a 
small modification is introduced. The fragment number of Ref. j3] was defined as 






(Al) 



where 




(A2) 



and 



fij F{dij , Qi) , fij Fyd^j, fl), fij F \dij , £ , fl) , 



(A3) 



dy = |Re(Z i -Z i )| 



(A4) 




(A5) 



g(k) = l + g e 



(k-M) 2 /: 



(A6) 



(A7) 
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In the present study, a gaseous nucleon is counted as a separate 'fragment' for the con- 
venience of the temperature measurement. We allow a few gaseous nucleons get spatially 
close by regarding the nucleon % as gaseous when the number of the neighboring nucleons 
including itself 

A 

is less than about three. This description is achieved by defining the fragment number Afp 
as 

A 

Ar F = J2^, J i = (l-w(ft))Xi 0) +w(q t ), (A9) 
i=i 

where 

w(q)=F(q,l/2a 2 q ,Q) (A10) 

and the parameters are taken as a q = 0.85 and Q = 2.5. The newly introduced function 
w(qi) becomes zero unless q < Q and the degree of isolation 2j reduces to the original one 
t\°\ On the other hand, Zj becomes close to one for the nucleons with qi < Q. 

The potential energies between the gaseous nucleons (<& < Q) are also modified in the 
present study. The term V mo d of Eq. (J3j) is introduced for this purpose. The potential energy 
between two gaseous nucleons is calculated by folding the widely spread wave functions with 
their interaction. The spread wave function, keeping the centroid position in coordinate and 
momentum space same as the original Gaussian wave function [Eq. ([!]) ]. is given by 



* a) = ilfiTJ ] exp 



3 



2u \ * 



(All) 

l + a v J v ; 

Z'i = ReZi +z(l + a)ImZi, (A12) 

where a is a smearing parameter. By denoting the potential energy between nucleon i and 
nucleon j as Vij calculated with the original Gaussian [Eq. ((TJ)] and Vij(a) calculated with 
the smeared Gaussian [Eq. (jAllj) ]. V mo d is given by 

A 

V mod = ^2w( y q i )w(q j ) -Vij+Vijia) . (A13) 

i<j 

The smearing parameter is chosen as a = 3 in order to reasonably reproduce the binding 
energies of small nuclei such as deuteron, triton and 3 He in case those stable compounds 
are formulated during dynamical processes. The density-dependent part of the effective 
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interaction is also included in V mo d with a similar approximation used in the triple-loop 
approximation [32]. When we evaluate Vy(a) with the density-dependent part, we have also 
modified the density p(r) consistently with the transformation of the wave function [Eq. 
JSDJ] as 

3 

p( r; a) = f (-Y e-^ a ^-^ 2 p(v')d 3 r'. (A14) 
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APPENDIX B: TEMPERATURE CALCULATION BY GASEOUS NUCLEONS 

The microcanonical temperature is defined by T _1 = (dS/dE)\y jA . We will show that 
this quantity can be evaluated with a good precision by using the average kinetic energy of 
a set of nucleons which can be chosen arbitrarily as long as two conditions are satisfied. The 
first condition is that the quantum effect should be negligible for any of these nucleons. The 
other condition is that these nucleons have to be chosen based on only the nucleon spatial 
coordinates without using momentum variables. In this section, we name the subsystem of 
these nucleons as System 1 and the rest of the system as System 2. 

The microcanonical ensemble {E, V, A} can be divided into two sub-ensembles: one with 
the number of nucleons of System 1 (Ax) is larger than zero and the other with A\ = 0. The 
density of microstates of the total ensemble is given by the sum of the microstates of these 
subensembles; 

W(E) = W Al>0 (E) + W Al=0 (E). (Bl) 

The constraints on the total volume V and the total nucleon number A are omitted in this 
expression and the following for brevity. The temperature of this ensemble is defined by 

1 d\nW(E) 



T dE 

W Al>0 (E) 1 W Al=0 (E) 1 



(B2) 

W(E) T Al>0 ' W(E) T Al=0 ' (B3) 

where T Al>0 = (<9(ln W M>Q {E) /dE))' 1 and T Al=0 = (<9(ln W Al =0 (E) / dE) ) _1 are the mi- 
crocanonical temperatures of the subensembles of W Al>0 and W Al=0 , respectively. The 
contribution of the second term of Eq. (jB3|) may be neglected if W Al =o(E) is much less than 
W Al>0 (E). 

Let us denote the spatial coordinates of nucleons in System 1 by 

i?! = {ReZi/v^; i E System 1}. (B4) 
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If the quantum effects are negligible in System 1, the density of microstates of System 1 
under the constraint of the nucleon positions R± is given by 



where E\ is the energy of System 1 (excluding the interaction energy with System 2), U(R\) 
is the potential energy within System 1 and M is the nucleon mass. The number of nucleons 
Ai is implicitly understood by R\ in the left hand side and in the following. Under the given 
condition Ri of the nucleon positions of System 1, we consider the density of microstates of 
System 2, which is denoted by W2(E' 2 , Ri). The nucleon positions of System 2 are constrained 
by the condition that only the the nucleons of Ri belong to System 1 when the defined 
algorithm is applied to the total system. System 2 consists of A — A\ nucleons. The energy 
E' 2 includes the interaction energy between System 1 and System 2 in addition to the internal 
energy of System 2 so that the total energy is E = E\ + E' 2 . The explicit expression for 
W 2 (E' 2 ,Ri) is not necessary. Using the density of microstates of System 1 and System 2, 
Wa^o^E) can be given by 



where the integral for Ri also includes the summation over the various cases of the nucleon 
number A\ (> 0) of System 1. 

Let us consider the case that Wa 1= o(E) is negligible compared with Wa^o^E) and thus 
W(E) is approximately equal to Wa^^E). Then the microcanonical temperature T can 
be calculated by 




(B5) 




(B6) 



T T, 



1 




Ai>0 
1 



1 



d\nW Al> o(E) 
dE 

f f d\nWx{E - E' 2 , Rx) 



xWx{E-E' 2 , R 1 )W 2 (E' 2 , R x )dE' 2 dR 



i 




(B7) 



By inserting Eq. (|B5|) to Eq. (|B7[). we obtain 




(B8) 
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where K\ = E x —U(Ri) is the total kinetic energy of System 1. It is also possible to evaluate 
the same temperature by 

T^im , (B9) 



3\A 



{E,A!>0} 



if the nucleons in System 1 follow the Maxwell distribution, which is usually the case with 
a good precision. In fact, when the distribution of the kinetic energy of System 1 is charac- 
terized by the Boltzmann factor e~ Kl ^ T °, Eq. (jB8|) can be calculated as 



U x - 1 \ !{U X - l)Kl M ~\- K ^dK x 



K i / {e,a 1>0 } J Kf Al ^e- K ^ T °dK x 

(U x - l^t^nu, - 1; 



= To 1 . (BIO) 
On the other hand, Eq. ()B9j) can be calculated as 



2 /K x \ 2 jKl Al e~ K ^dK x 



3 \ A i I {e,Ai>o} 3Al f Kl M ~ X e~ K ^dK x 



3 



2 T 2 _ Y{\A X + 1) 

^ tI m y{\a x ) 

= T . (Bll) 

Equation (|B8|) is derived under the assumption that Wa 1= o(E) fW{E) is negligible. If 
Wa 1= q(E) is not negligible, the applicability of Eq. (|B8J) and Eq. ()B9jl is questioned unless 
X^ 1>0 ~ Ta 1= o is satisfied. If we can not evaluate T^ 1=0 practically, Eq. (jB8|) and Eq. (jB9|) 
should be used carefully when the average number of the nucleons in System 1 is less than 
one. 
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